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Abstract 

We solve some forms of non homogeneous differential equations in one 
and two dimensions. By expanding the solution into whell-posed closed 
form-Eisenstein series the solution itself is quite simple and elementary. 
Also we consider Fourier series solutions of linear differential operator 
equations. In the third section we study operators which are functions of 
the Leibnitz derivative. The last result is the complete solution of a non 
homogenus 2-degree ODE with linear coeficients. The non homogenous 
part is an arbirtary function of Z/2(R). 



1 The Divisor Sums and ODE 

Proposition 1. If x, is positive real number and / is analytic in (-1,1), 
with /(O) = 0, then 

/fx \ oo f( d Uo) 

exp( / /(Odt) =Y[(l- e -n^-*-" n/d) (1) 

Woo / n=1 

Where n is the Moebius function. See also and [Ap] chapter 2 for the 
Moebius as also for other multiplicative functions. 
Proof. See [B] 



Examples of such identities are 

1) Let = l[n] = 1 if n = 1, else. X(n) = -fx(n) 



new 



»«(») 
|— =e -« 



2) Let 1^21 = „, n = 1,2,... then f(x) = and 



1 



Where 4>{n) is the Euler Totient function. 

Proposition 2. If A(n) is arbitrary sequence of numbers we have for 
x > 

dT (f, T^n Md)Kn/d) \ ~ Edm A(d)(-dY^n/d) 
dx v \ 2-~t e nx - 1 / e nx - 1 

\n = l / 71=1 

Proof. 

See also [B] 

We will use Proposition 1 to find the solution of the N — th degree linear 
differential equation 

^ <r - E dk g(dV(n/d) - 



Lemma 1. Set 

JV 

P(i) = ^a„i" (4) 

i/=0 

then the solution of (3) is 

n — 1 n — 1 

Proof 1. From Proposition 2, it is clear that if 



e" 1 - 1 



n = l 

for a certain A(fc), then (3) becomes 

AT /no v-»oc , / ,x i i ,n \ (") 



V a M f V ^MrVrf) y _ y> Ed|n g(rf)M(n/rf) 

i/=0 \n = l / n = l 



or 



^ , E„|„ A(d)(-d)V(n/d) ~ E„| B C(dUn/d) 

} CLvix) > = > 

/ j v ' / j gni ^ / ^ gnx ^ 

f — n— 1 rt — 1 

or 

E d |„^)-P(-rf)Mn/rf) _ - E d |„^)M»/d) 

n — 1 n — 1 



Hence it must be 



A(d)P(-d)fi(n/d) = Y C(d)fi(n/d) (6) 

d\n d\n 
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Eq.6 shows clearly that 

A(n) = C(n)/P(-n) 
Also one can see that we have 



W(n) 

qTLX ' 

n — 1 n= 1 \ d|n 

and 



2_. = 2^ ^< n ) e 

n— 1 n— 1 

Proof 2. Let 

oc 

/ \ \ ^ — rix 

u(x) = u n e 

n=l 

then setting into (3) the above expansion we get the same result in a more 
easy way. 

Theorem 1. If 

N 00 n( \ 

y av u^u{ X ) = y-^- 

Then 



P*(n) 

n — i 

, where 

JV 
fe=0 

Proof. As in Lemma 1. 

Next we proceed with the 2-dimension problem with a similar way. 



We set 

N 



fc=0 

and 

M 



k=0 

M 

P y (w) = J2biw l (8) 



Also 

oo 
k,m — 1 
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Theorem 2. The equation 

N,M 

^a k b l u W ' (l \x,y) = G(x,y) (10) 



k,l=0 

have solution 



\— S(n,m) 

u ^»)= E (en,.!)^,!) ( U ) 

S(n,m) = ^2 B(d,S)fi{n/d)fj,(m/8) (12) 

d\n,S\m 

B ^ m)= P x {-n)P y {-m) E C ^ < 13 > 



n,m=l 

where 



d|n,(5|m 

and 



fe|n,r|m 

Proof. Let 

, v 2 ^ Sdin A d (y)n(n/d) 
u(x, y) = ^ ^73^ 

n=l 

Then differentiating with respect to x we get 

,<'»Vv)=x; E ^-f ( " /d) w 

n=l 

then with respect to y we get 

n=l 

but 

A *(v) = 2^ (c) 

combining the above we get the result. 



Note. The polynomials P apparently must have no solutions in natural 
numbers. 



Examples. 

1) Set 



P x (w) = 1 + V2w + w 2 

also set 



Py(w) = 1 + V2W 

Then the equation 

u(x,y) + V2u x (x, y) + u xx (x,y) + V2u y (x, y) + 2u xy (x,y) + V2u xxy (x, y) 
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oo oo 

' / e na: — 1 ' ' e my — 1 



have solution 



where 



2) If 



V U> \ ^ 1 - V2n + n 2 \ ^ 1-V2r, 

\n= 1 / \m— 1 



d\r. 



P x (w) = 2 + V3w + ', 



2 

- W 



Py(w) = 1 + VSW + ' 

then the solution of 

2li ~f" \/^Ux ~\~ Uxxx ~~\~~ 2\/3liy ~\~ t^Uxy ~\~ ^SUxxxy ~t~ 2lZyy ~\~ "*J~^>Uxyy ~~\~~ Ux 

no 

log(n + m) 



E 



(e nx - l)(e m y - 1) 



u{x,y)- 2^ ( e ™ - l)(e™v - 1) 

n,m= 1 

E d |„E 5 | m 1 °g( d + ' 5 ) 



(c" 1 - l)(e m » - 1 

n,Tn=L 

where 

B(n, m) 

(2-^n-n 3 )(l-\/3m + m 2 ) 

Observe that in this example we are not able to split the solution into two 
parts in x and y. 

Note. For no confusion the form of the equation is defined by 

N,M 



J2 a n b m u^' {m \x,y) 



n,m=0 

where the a„ and b n are respectively that of P x and P y . 

2 Series Solutions 

If F is an operator such that 

(t\) k d k 



k\ dx k 

k=0 



We will try to solve the equation 
d 



dt u(x,t) = -u{x,t) 
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Assume that 

u(x,t) = fix) 



then 



= e«*^/(x) = ^e**/(») = ^«(x >t ) 



Hence the operator exp {t-§^) produces the solution. From Eq.(3) and 
Lemma one can take the limit N — > oo then 

oo 

F (^) u(x, t) = exp (-t-^j u(x, t) = f{x) = J2 C(n)e~ nx 



thus according to Lemma the solution of the above equation must be 



Using the parameter A which can take and complex values one can arrive 
to the conclusion that 

J^7 M (M) + J^«(M) = (14) 

have solution 

u(x,t) = Y ^ -e~ nx (15) 

n— 1 

One can see that, a solution of Schrodingers equation 

Qm 

«t(M) = -^«(a;,«) + V(t)u(a;,«) (16) 



«(M) = > TT-^T 777^7 ( 17 

Z-j exp (-n)"*t - /(t) 

n— 1 

where 

f{x) = Viw)dw (18) 
This is the general case in which the potential depends only in time. 



Let again 

JV 

P(i):=^o k i' (19) 



fc=0 

Consider now the equation 



N 

E w(t) w = p (3* )=I_I (20) 



G 



the solution is 



y(x) = Xao - a " 2 n - ai + V C{k)e"» (21) 
a ° k=i 

Where p^ is the roots of P(x) = The same equation have solution 
according to the Theorems of section 1: 

y ( x ) = * e <"» (22) 

yy > A^i nPUn) v ' 

nez* 



An interesting question is how one can extract from (21) and (22) the 
roots pk- 

Anyway when if we let N — > oo, then 

Theorem 3. 





nF(in) 
nez* 



ii) 
iii) 



F (^) y(x) = x - 7T (24) 



y(x) = ~ "'^ ' + V- C(fc)e p^ (25 ) 
a ° fe =i 

the pfc are roots of Eq. F(x) — x — it. The function F must have not 
integer roots in the imaginery line. 



Now consider the function F(x) = e 27r%x +x — ir — l. It is a = F(0) = -7r 
and ai = — 2m + 1 and pfe = k. Hence the two representations are 



iez" 



n(e 2nn +in-n-l) 



I n -a;7T + 7T +27TI-1 >r-^ 



n— — oo 



the differential equation is 

y(iE — 2m) + y'(x) — y(x)(-K + 1) = x — 7T 

Examples 
1) 

F(z) = + a; + 1 

then 

VW= E w (-n + l+ - C 
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and also 

oo 

/ \ 1 71 ~\- X ttix \ ^ s~i / 1 \ 2k7vix 

y{x) = g +e 2^ c ^ e 

fc — — oo 

in this example the diferential eq. is 

y(x + 1) + y'{x) + y(x) = x - n 

2) For F(x) = cosh(;r) + x + 1 the DE is 

y{x + 1) + y(x - 1) , 

- hy(i)+ t/(a;) = x - tv 

with solution 



The series in a first view can not become more fast convergent. If we 
consider for example 

F(x) = cos(a;) + x + 1 

then 

y M (x)= n ( cosh (n) + m+l) eln3: 

|n|<M V W ' 

But for the diferential equation holds 
y M {x + + y M (x - i 



+ v'm(x) + y M (x) -x + 7: = 



2 

If we try with F(x) — cos(nx) + x + 1 then we have 

i 



VM 



(*)= E 



n(cosh(n7r) + imr + 1) 

|<m x x > > 



and for the diferential equation holds 

y M {x + in) + y M {x-m) , (I 

2 + Vm( X ) +yM{x)-X + TT = 0^ — 

which is the same. 

Thus we can say that even the equations are not easy to solve numericaly, 
the solutions itself under certain conditions may behave very good i.c 
ynf(x) is very fast convergent. 

3) Let us consider now a curius case. The I/2(Ft) function h(x) = e x ~ e 
then if F(x) = h(x) + x + 1 (we dont need the roots), we have 

" (-1)' 



E /T -y(x + l + l) + y'(x) + y(x) = x - n (26) 
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and 



y(x) = V — — e inx (27) 

^ n(m + 1 + e™-« 1 ") 

is the solution indeed. 
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3 Functions of ^ 

We proceed with the following 
Lemma 2. Let 

oc 
n=0 

be analytic function in R such that for every a, b > there exist constant 
depending from f,Mf. 

\f(x)\<M f (l + \xne-"W (28) 

Let also <f>(x) real valued function with values in C, such that for every 
c > there exist costant M$: 

\4>{x)\ <M |x| c (29) 

then 

Jo fc=0 fc=0 

(30) 

Proof. 

See [Ba] 

If happens L(f> — y(x) then 

f-OO 

/ f(x)(L ( - 1) y)(x)e- xa dx = 
Jo 

_ /( 2fc > (0) d^y(s) /P*+^> (0) rf 2fc +y g ) 

^ (2fe)! ds 2fc ^ (2fe + l)! ds 2 ^ 1 1 j 

fc— fc— 

hence we can write 

r f(x)(L^y)(x)e- xs dx = f e (^) »( a ) - /„ (^) y( S ) 

«/ 

f(-x)(L^y)(x)e- xs dx = / (^) y(s) (32) 

Theorem 4. 

i) It holds 

/ (£) w ( a ) = L (/(-x)^- 1 ^) (*)) (a) (33) 

ii) The solution of 

/ {£) Vis) = 9{s) (34) 

is 

y { s) = j~^Me--dx (35) 
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This theorem shows clearly that we can find solutions in integral-closed- 
form, of the general not homogeneous equation (if existing the Laplace 
transforms) . 



Examples. 

1) If hapens g(x) = 1/x 2 , then (L ( -~ 1 ^g)(x) = x and thus the solution of 
/ (|) y« = l/, 2 (36) 



y( x ) = J j^x) e ' xadx ^ 

where the form of the equations (36) and (37) is that of (34) and (35). In 
the same way as in the above example we can set other values for g(x) 

2) Set 

OO 

h(x) = ^c k x Vk e- lkX 

k=0 

Then the equation 

oc 

Y,c k y M (x-l k )=g(x) (38) 

k=0 

have solution 

= / ^-T7 — T-^e dw 39 
Jo h (-w) 

This method is like solving (38) with Fourier or Laplace transforms but 
we avoid some restrictions of y and g to be in Z,2(R). Note also that it is 
solved with Laplace theory. 

3) We try now to evaluate T = 1 _^j L _ ■ Let Ty(x) = g(x), then 
1 l d_ + jf_ (f | 



1 + da; da; 2 dx 3 

dx 



1 „ _ i d ^ , d2 v dS y 



^= T y^ = YTj y{x) = 1 -fx + ^-^ + 

dx 

we use (39) and get 

/■oo 

j/(a;) = / (Z/ _1 'g)(w)(l — w)e~ xw dw 
Jo 

solving with respect to g we have 

g(w) = Ty(w) = [°° e~™dx 



l-x 



In general holds 



h {i) f{x) = j (L^ftMH-vJe'^dw (40) 
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under certain conditions of convergence. For example if y(x) =polynomial 
in x. Also for h(x) = log(l + x), then 



n=l 

/■ OO 

= / log(l — w)e~ xw dw 

Jo 

Thus for example if we consider the equation 

f = (e) 

with /(x) = x + log(l + x), then the differential equation (e) is actualy 
the 



[ (L^- 1) y)(w)log(l-w)e- xw dw = g(x) 
dx Jo 



and have solution 



J Q -w + log(l - w) 



Relation (40) is very useful if one can set a one to one relation between 
h and a function of y. For example if one take h(x) — e x then for all the 
functions 

Q(x) + e x 

with Q(x) = X^fcLo akxk > tne equation will be 

N d k 

^2 ^dx*^^ + y( - x + ^ = 9 ^ 

k=0 

will have the same solution type 



Can we say that for a given function /i(a;), (e 1 in the examples) exists a 
unique form of y(x) (such as y(x + 1) in the examples)? The exponential 
and polynomial functions behave very good, but what happens with other 
values of h. For example: Exist h 2 (x) giving us the form y(x) 2 in the 
differential equation? Then all the differential equations, with this term, 
will be solvable in integral forms with knowing one function only. But 
this seems not to happen. It happens with the differential functions, i.e. 

fc(£)v(*) = iog(i-£). 

For example if we consider the equation 

2 h y{x) " a -ah v{x) - log ( : 1) y{x) = 9{x) 
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the solution is 



, / ( L( 1] 9)( w ) e -™ dw 



Another related equation is 

d 



which have a solution 



y(x + + log ^1 - — ^ y(x) = g(x) 

f 



' e~' w + log(l + w) 



From the above examples and (39) and (40) one can see that the inversion 
with respect to some g(x) is 
Theorem 5. (Inversion) 



l U) h (£) 



(41) 



v dx 

This means if 



y(x)= [ y ~ u( y /^' e- xw dw (1.2) 



ft(-io) 
then 

/>oo 

y ( - 1) (x)=g(x) = (L { - 1) y)(w)h(-w)e- xw dw (43) 
Jo 

Where 



&(^)y(aO = s(aO = S/i(aO 



and 

Example. 

If h(x) = e~ x + 1 then fr'" 1 '^) = -L = -r^. This means that if 
g(x) = l/x 2 

y(x + 1) + y(x) = g(x) 

then 

K .,.«(_*( u + :) + *( 1 ,*±l)) 

and 

poo 

g(x) = yi(x)= / (i (_1) 2y)(w)(e"' + l)e~ xw dw = 1/a; 2 



Where ^ is the Polygamma function i.e 

i dr( 2 ) 

and 

V»(n, 2 ) = V (n) W 
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(see and Mathematica notes). The above example is trivial and can be 
solved with Laplace theory. 



Now we will find a way to solve the equation 

[aix + bi)f"(x) + (a 2 x + b 2 )f'(x) + (a 3 x + b 3 )f(x) = g{x) (44) 
where /, g e L 2 (R). 

Let the Fourier Transform of a function of L 2 (R) is 

/oc 
f(t)e-^dx 
-CO 

the Inverse Fourier Transform is 

/(*) = ^ r tw*^ 

J — 00 

Lemma 3. 

/OO 
f{x)x n e- ix ' 1 dx = i n (f) {n \ 1 ). (45) 
- OO 

(/W)( 7 ) = (i 7 )"7(7). (46) 

/OO 
/'(^(a;)^ 7 ^ = 
-OO 

/oo /*oo 
/(^'(aOe^dz + H-y) / /(a:)A(x)e~ la;7 da;. (47) 
-co J — oo 

/oo 1*00 
/"(rc^Or)^ 7 = / f(x) A" (x)e~ lxl dx+ 
- oo J — oo 

p co poo 

2(-»7) / /(aOv4>)e _i * 7 da; + (-z 7 ) 2 / /(^(aOe"^ 7 ^. (48) 



Proof. 

The proof of (45) and (46) can be found in [Pa]. The relations (47) and 
(48) are obtained with integration by parts. 
Theorem 6. 

When f,gE I/2(R) and lim^^^ |/(a;)a; 2+t | = 0, e > 0, equation (44) 
can reduced in into 

f l \ 

(-ia 1 j 2 +a2j+ia 3 )— h(-fei7 2 -2iai 7 +i62 7 + a 2+&3)/(7) =d{l) (49) 

a 7 

which is solvable. 
Proof. 

Take the Fourier Transform in both sides of (44) and use Lemma 3. 
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